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Abstract 

We develop a scheme to construct the Hamiltonians of the lambda, vee and 
cascade type of three-level configurations using the generators of SU (3) group. It 
turns out that this approach provides a well defined selection rule to give different 
Hamitonians for each configurations. The lambda and vee type configurations are 
exactly solved with different initial conditions while taking the two-mode classical 
and quantized fields . For the classical field, it is shown that the Rabi oscillation 
of the lambda model is similar to that of the vee model and the dynamics of the 
vee model can be recovered from lambda model and vice versa simply by inversion. 
We then proceed to solve the quantized version of both models introducing a novel 
Euler matrix formalism. It is shown that this dynamical symmetry exhibited in 
the Rabi oscillation of two configurations for the semiclassical models is completely 
destroyed on quantization of the field modes. The symmetry can be restored within 
the quantized models when the field modes are both in the coherent states with 
large average photon number which is depicted through the collapse and revival of 
the Rabi oscillations. 
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I. Introduction 



Quantum Optics gives birth to many novel proposals which are within reach of present- 
day ingenious experiments performed with intense narrow-band tunable laser and high-Q 
superconducting cavity [1]. Major thrust in the atomic, molecular and optical experiments 
primarily involves the coherent manipulation of the quantum states which may be useful to 
verify several interesting results of quantum information theory and also the experimental 
realization of the quantum computer [2,3]. The actual number of the quantum mechanical 
states of atoms involved in the interaction with light is of much importance in these days 
since many coherent effects are due to the level structure of the atom. It is well-known 
that the two-level system and its quantized version, namely, the Jaynes-Cummings model 
(JCM), have been proved to be an useful theoretical laboratory to understand many 
subtle issues of the cavity electrodynamics [4,5]. The two- level system is modeled using 
the Pauli's spin matrices - the spin-half representation of SU (2) group, where apart from 
the level number, the spectrum is designated by the photon number as the quantum 
number. A natural but non-trivial extension of the JCM is the three-level system and 
it exhibits plethora of optical phenomena such as, two-photon coherence [6], resonance 
Raman scattering [7], double resonance process [8], population trapping [9], three- level 
super radiance [10], three-level echoes [11], STIRAP [12], quantum jump [13], quantum 
zeno effect [14], Electromagnetically Induced Transparency [15,16] etc. There are three 
distinct schemes of the three-level configurations which are classified as the lambda, vee 
and cascade systems respectively. The Hamiltonians of these configurations are generally 
modeled by two two-level systems coupled by the two modes of cavity fields of different 
frequencies [17,18]. Although these Hamiltonians succeed in revealing several phenomena 
[19,20], however, their ad hoc construction subsides the underlying symmetry and its 
role in the population dynamics of these systems. The connection between the SU(N) 
symmetry and the iV-level system in general, was investigated extensively in recent past 
[21-27]. These studies not only mimic the possible connection between quantum optics 
with the octet symmetry, well-known paradigm of particle physics, but for N — 3, it also 
reveals several interesting results such as the realization of the eight dimensional Bloch 
equation, existence of non-linear constants [18,22], population transfer via continuum 
[28], dynamical aspects of three-level system in absence of dissipation [29] etc. However, 
inspite of these progress, a general formalism as well as the ab initio solutions of all three 
configurations are yet to be developed for the reasons mentioned below. 

The generic model Hamiltonian of a three-level configuration with three well-defined 
energy levels can be represented by the hermitian matrix 
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where hij(i,j = 1, 2, 3) be the matrix element of specific transition and Aj be the detuning 
which vanishes at resonance. We note that from Eq.(l), the lambda system, which corre- 
sponds to the transition 1 <-> 3 <-> 2 shown in Fig. la, can be described by the Hamiltonian 
with elements h 2 ± = 0, h 32 ^ and h 3 i ^ 0. Similarly the vee model, characterized by 
the transition 3 <-> 1 <-> 2 shown in Fig. lb, corresponds to the elements h 2 \ ^ 0, /i32 = 
and h 3 \ ^ and for the cascade model we have transition 1 <-> 2 <-> 3, we have /121 7^ 0, 
^32 7^ and /i3i = respectively. Thus we have distinct Hamiltonian for three different 
configurations which can be read off from Eq.(l) shown in Fig.l. This definition, how- 
ever, differs from the proposal advocated by Hioe and Eberly, who argued the order of 
the energy levels to be E\ < E 3 < E 2 for the lambda system, E 2 < E 3 < E\ for the 
vee system and E\ < E 2 < E 3 for the cascade system respectively [18,21,22]. In their 
scheme, the level-2 is always be the intermediary level which becomes the upper, lower 
and middle level to generate the lambda, vee and cascade configurations respectively. It 
is worth noting that, if we follow their scheme, these energy conditions map all three 
three-level configurations to a unique cascade Hamiltonian described by the matrix with 
elements h± 2 7^ 0, h 23 7^ and h 13 = in Eq.(l). Thus because of the similar structure 
of the model Hamiltonian, if we start formulating the solutions of the lambda, vee and 
cascade configurations, then it would led to same spectral feature. Furthermore, due to 
the same reason, the eight dimensional Bloch equation always remains same for all three 
models [18,22]. Both of these consequences go against the usual notion because wide 
range of coherent phenomena mentioned above arises essentially due to different class of 
the three-level configurations. Thus it is worth pursing to formulate a comprehensive 
approach, where we have distinct Hamiltonian for three configurations without altering 
the second level for each model. 

The problem of preparing multilevel atoms using one or more laser pulses is of con- 
siderable importance from experimental point of view. Thus the completeness of the 
study of the three-level systems requires the exact solution of these models to find the 
probability amplitudes of all levels, the effect of the field quantization on the population 
oscillation and, most importantly, the observation of the collapse and revival effect. In 
recent past, the three- level systems and its several ramifications were extensively covered 
in a general framework of the SU(N) group having iV-levels [21-27,30,31]. Also, the semi- 
classical model [24,32,33] and its fully quantized version [23,34,35] are studied, but to our 
knowledge, the pursuit of the exact solutions of different three-level systems in the spirit 
of the theory of Electron Spin Resonance (ESR) model and JCM, are still to be facilitated 
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analytically. 

In a recent paper, we have studied the exact solutions of the equidistant cascade sys- 
tem interacting with the single mode classical and quantized field with different initial 
conditions [36]. It is shown that for the semiclassical model the Rabi oscillation exhibits 
a symmetric pattern of evolution, which is destroyed on quantization of the cavity field. 
We also show that this symmetry is restored by taking the cavity mode to be the coherent 
state indicating the proximity of the coherent state to the classical field. We have further 
studied the equidistant cascade four- level system and obtain similar conclusions [37]. To 
extend above studies for the lambda and vee models we note that the vee configuration 
can be obtained from the lambda configuration simply by inversion. However, it is worth 
noting that, the lambda configuration is associated with processes such as STIRAP [12], 
EIT [15,16] etc, while the vee configuration corresponds to the phenomena such as quan- 
tum jump [13], quantum zeno effect [14], quantum beat [3] etc indicating that both the 
processes are fundamentally different. It is therefore natural to examine the inversion 
symmetry between the models by comparing their Rabi oscillations and study the effect 
of the field quantization on that symmetry. The comparison shows that the inversion 
symmetry exhibited by the semiclassical models is completely spoiled on quantization of 
the cavity modes indicating the non-trivial role of the vacuum fluctuation in the symmetry 
breaking. 

The remaining part of the paper is organized as follows. In Section-II, we discuss 
the basic tenets of the SU(3) group necessary to develop the Hamiltonian of all possible 
three-level configurations. Section-Ill deals with the solution of the lambda model taking 
the two field modes as the classical fields and then in Section IV we proceed to solve 
the corresponding quantized version of the model using a novel Euler matrix formalism. 
Section-V and VI we present similar calculation for the vee model taking the mode fields 
to be first classical and then quantized respectively. In Section- VII we compare the 
population dynamics in both models and discuss its implications. Finally in Section- VIII 
we conclude our results. 

II.The Models 

The most general Hamiltonian of a typical three-level configuration is given by Eq.(l) 
which contains several non-zero matrix elements showing all possible allowed transitions. 
To show how the SU(3) symmetry group provides a definite scheme of selection rule which 
forbids any one of the three transitions to give the Hamiltonian of a specific model, let us 
briefly recall the tenets of SU (3) group described by the Gell-Mann matrices, namely, 
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These matrices follow the following commutation and anti-commutation relations 

[Ai, \j] = 2ifijk\k, {Aj, \j} = + 2dijk\k, 



(3) 



respectively, where and = 1,2,. .8) represent completely symmetric and 

completely antisymmetric structure constants which characterizes SU(3) group [39]. It is 
customary to define the shift operators T, U and V spin as 



T± = i(Ai±iA 2 ), C/± = |(A 6 ±iA 7 ), 
They satisfy the closed algebra 



V± = !(A 4 ±iA 5 ). 



[U+,U-] = U 3 , 




= v 3 , 


[T + ,T_]=T 3 , 


[T 3 ,T ± ] = ±2T ± , 


[T 3 ,U ± 


= tu±, 


[T 3 ,V±] = ±V±, 


[V 3 ,T ± ]=±T ± , 


[v 3 ,u ± ] 


= ±u ± , 


[V 3 ,V ± ] = ±2V ± , 


[U 3 ,T ± } = T T ± , 


[u 3 ,u ± ] 


= ±2U±, 


[U 3 ,V ± ] = ±V ± , 


[T + ,V-] = -£/_, 


[T+, U + 


] = v + , 


[U+,V-]=T., 


[T_,V + ]=U + , 


F-,u.]-- 


= -V-, 


p.,V + ] = -T+, 


where the diagonal terms 


are T 3 = 


A 3 , U 3 = 


(V3\ 8 - A 3 )/2 and V 3 



(4) 



(5) 



respectively. 

The Hamiltonian of the semiclassical lambda model is given by 

H A = Hf + Hf r , (6a) 
where the unperturbed and interaction parts including the detuning terms are given by 
Hf = uo 2 )V 3 + h(n 2 -lui- ^ 2 )T 3 , (6b) 

and 

H£ = h(AfV 3 + A^T 3 )+ 
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hKi(V + exp(— iQit) + V- exp(iQit)) + Tik 2 {T + exp(— iQ 2 t) + T- exp(iQ 2 t)), (6c) 

respectively. In Eq.(6), Q t {% = 1, 2) are the external frequencies of the bi-chromatic field, 
Ki are the coupling parameters and hu>i(— — Ei),tno 2 (= —E 2 ),h(u 2 + ^i)(— E 3 ) be the 
respective energies of the three levels. A^ = (2uo x + uo 2 — Q\) and A2 = (0J1 + 2oo 2 — Q 2 ) 
represent the respective detuning from the bi-chromatic external frequencies as shown in 
Fig.l. 

Proceeding in the same way, the semiclassical vee type three-level system can be 
written as 

H v = H7 + Hg, (7a) 

where 

h7 = h(Q! - u! - u 2 )v 3 + %{Q 2 ^ 2 )u 3 , (7b) 

and 

%Ki(V + exp(— iQ.it) + V- exp(iQit)) + %k 2 (U + exp(— iQ 2 t) + U- exp(iQ 2 t)) (7c) 

where A] 7 = (2coi +uo 2 — Qi) and A^ = (2oo 2 + uoi — Q 2 ) be the detuning shown in Fig.2. 
Similarly the semiclassical cascade three-level model is given by 

H s = Hf + fig, (8a) 

where 

Hf = h{Qi +u 2 - ui)\J 3 + h{Q 2 + cui- lu 2 )T 3 , (8b) 

and 

Hf I = ft(Af[/3 + AiT 3 )+ 

%Ki(U + exp(-iQit) + U- exp(iQit)) + hn 2 (T + exp(-iQ 2 t) + T_ exp(iQ 2 t)) (8c) 

respectively with respective detuning Af = (2ui — u 2 — Qi) and A2 = (2u 2 — u>i — Q 2 ). 

Taking the fields to be the quantized cavity fields, in the rotating wave approximation, 
the Hamiltonian of the quantized lambda configuration is given by 

H A = + Hfj, (9a) 

where, 
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Hf = h(tt 2 - wi - uj 2 )T 3 + hiQi - U! - co 2 )V 3 + E fyaUj, 



Hjj = hAfV 3 + HA% T 3 + ^(F+ai + V_a\) + % 2 (T+a 2 + T_a|), 



(9b) 



(9c) 



where af and a, (i = 1, 2) be the creation and annihilation operators of the cavity modes, 
gi be the coupling constants and fli be the mode frequencies. Proceeding in the similar 
pattern, the Hamiltonian of the quantized vee system is given by 



tf 1 



H? + HY Z , 



where, 



Hj v = h(n 2 -uj 1 - lu 2 )XJ 3 + - wi - co 2 )V 3 + e fy<4 

3=1 



HjY = hA\V 3 + hA^V 3 + h gi (V +ai + V-a{) + hg 2 (U + a 2 + 

respectively. Similarly the Hamiltonian of the quantized cascade system reads 
H~ = Hf + Hfj, 

where 

2 

#f = ft(fi 2 — o;i — w 2 )T 3 + K(£li —ui— co 2 )U 3 + Qja)aj, 

3=1 

H§ = HAfU 3 + ftA|T 3 + h gi (U +ai + U-a[) + hg 2 (T + a 2 + T_a\). 



(10a) 

(10b) 
(10c) 

(11a) 

(lib) 
(11c) 



Using the algebra given in Eq.(5) and that of field operators, it is easy to check that 



[Hf, H l n } = for Aj = -Ai 



A and V) for the lambda and vee model and A^ = A2 



for the cascade model which are identified as the two photon resonance condition and 
equal detuning conditions, respectively [18,21,22,24,26]. This ensures that each piece 
of the Hamiltonian has the simultaneous eigen functions. Thus we note that, unlike 
Ref. [18,21,22], precise formulation of the aforementioned three-level configurations require 
the use of a subset of Gell-Mann Aj matrices rather than the use of all matrices. We now 
proceed to solve the lambda and vee configurations for the classical and the quantized 
field separately. 



III. The semiclassical lambda system 
At zero detuning the Hamiltonian of the lambda type three-level system is given by 

h(ui + uj 2 ) hn 2 exp[—ifl 2 t] %K\ exp[— iVL\t] 
II ' ///, L .cxp HI, I -hu 2 . (12) 

Hki exp[iQit] — Hcui 
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The solution of the Schrodinger equation corresponding to Hamiltonian (12) is given by 

= d(t) |1) + C 2 (t) |2) + C 3 {t) |3) (13) 

where C\(t), C 2 (t) and C 3 (t) be the time-dependent normalized amplitudes of the lower, 
middle and upper levels with the respective basis states, 





1 



|2> 




1 




|3> 



1 






(14) 



respectively. We now proceed to calculate the probability amplitudes of the three states. 
Substituting Eq.(13) in Schrodinger equation and equating the coefficients of |2) , |3) and 
|1) from both sides we obtain 



' at 



(uo 2 + oo\)C 3 + K X exp(-ifiit)Ci + k 2 exp(—iQ 2 t)C 2 , 



i^jf = —oj 2 C 2 + k 2 exp(i£l 2 t)C 3 , 
= -uiCi + Ki exp{iVlit)C 3 . 

Let the solutions of Eqs.(15a-c) are of the following form, 
Ci = A 1 exp(iS 1 t), 
C 2 = A 2 exp(iS 2 t), 
C 3 = A 3 exp(iS 3 t), 



(15a) 
(15b) 
(15c) 

(16a) 
(16b) 
(16c) 



where A^s' are the time independent constants to be determined. Putting Eqs.(16a-c) in 
Eqs.(15a-c) we obtain 



(S 3 + uj 2 + cji)A 3 + k 2 A 2 + /tiAi = 0, 
(S 3 + tt 2 - uj 2 )A 2 + k 2 A 3 = 0, 
(S 3 + fii - wi)Ai + /tiA 3 = 0. 



(17a) 
(17b) 
(17c) 



In deriving Eqs.(17), the time independence of the amplitudes A 3 , A 2 and Ai are ensured 
by invoking the conditions S 2 = S 3 + Q 2 and Si = S3 + Qi. At resonance, we have 
= = -A£ i.e, (2u 2 + u ± ) - Q 2 = = (u 2 + 2wi) - fti and the solution of Eq.(17) 
yields 



S 3 = -(u 2 + u 1 )±A, 



(18a) 
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S 3 = -(w 2 + wi) (18b) 
where A = \fK%~+~K% and we have three values of 52 and Si namely 

Si = tu 2 , S 2 ' 3 = u 2 ± A, (19a) 

S-i = Wl> Si' 3 = U! ± A. (19b) 

Using Eqs.(18) and (19), Eq.(16) can be written as 
C 3 (t) = Al exp(-i(cu 2 + u ± )t) 

+A|exp(i(-(cj 2 + cj 1 ) +A)t) +A|exp(i(-(cu 2 + cJi) - A)t), (20a) 

C 2 (t) = Al exp(iu 2 t) + A 2 2 exp(i(iu 2 + A)t) + Af (i(w 2 - A)t), (20b) 

d(t) = A} exp(iwit) + A 2 exp(i(cjx + A)t) + A?(i(wi - A)t), (20c) 

where A r s are the constants which can be calculated from the following initial conditions: 
Case-I: At t = let the atom is in level-1, i.e. d(0) = 1, C 2 (0) = 0, C 3 (0) = 0. Using 
Eqns (15) and (20), the corresponding time-dependent probabilities of the three levels are 

\C 3 (t)\ 2 = ^ sin 2 At, (21a) 

|C 2 (t)| 2 = 44^sin 4 At/2, (21b) 

|Ci(t)| 2 = ^t(«l + /c? cos At) 2 . (21c) 

Case-II: If the atom is initially in level-2, i.e. d(0) = 0, C 2 (0) = 1 and C 3 (0) = 0, the 
probabilities of the three states are 

|C 3 (t)| 2 = ^sin 2 At, (22a) 

|C 2 (t)| 2 = ^(/t 2 + ^cosAt) 2 , (22b) 

|Ci(t)| 2 = 4^ sin 4 At/2. (22c) 

Case-Ill: When the atom is initially in level-3, i.e. d(0) = 0, C 2 (0) = and C 3 (0) = 1, 
the time evolution of the probabilities of the three states are 

|C 3 (t)| 2 = cos 2 At, (23a) 
|C 2 (t)| 2 = ^ sin 2 At, (23b) 
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|Ci(f)| 2 = $siii 2 Af. 
We now proceed to solve the quantized version of the above model. 



(23c) 



IV. The quantized lambda system 

We now consider the three-level lambda system interacting with a bi-chromatic quan- 
tized fields described by the Hamiltonian Eq.(9). At zero detuning the solution of the 
Hamiltonian is given by 



I*a(*)>= E PI 

n,m=0 



n— l,m+l 



(t) \n - 1, m + 1, 1) + CZ m (t) \n, m, 2) + C^ L ' m (t) \n - 1, m, 3)], 



(24) 



where n and m represent the photon number corresponding to two modes of the bi- 
chromatic fields. This interaction Hamiltonian that couples the atom-field states \n — 1, m, 3), 
| ra, m, 2) and \n — 1, m + 1, 1) and forms the lambda configuration shown in Fig.l is given 
by 



Hjj = h 



92\fn gi\/m+ 1 
g 2 Vn 
gi\ / m+ 1 



The eigenvalues of the Hamiltonian are given by \± = ±fry ngl + (m + l)g 2 
and A = 0(= f2o), respectively with the corresponding dressed eigenstates 



\nm, 3) 
| Tim, 2) 
|nm, 1) 



T n>m (gi,g 2 ) 



\n 



n — 1, m, 3) 

|n, to, 2) 
- l,m+l,l) 



(25) 
(26) 



In Eq.(26), the dressed states are constructed by rotating the bare states with the Euler 
matrix given by 



T„, m (gi,g2) 



C 3 C 2 - C1S2S3 C3S2 - C1C2S3 S 3 S 1 
~ S 3 C 2 — C1S2C3 — S3S2 + c lC 2 C3 C3S1 
S\S 2 ~ s l c 2 C\ 



(27) 



where Sj = sin#j and q = cos#j (i = 1, 2, 3). The elements of the matrix are found to 

71 92 V 2(ngf + (m+l)gf ) 5,1 V 2(n S 2+( m+ l) 9 2) 

T n>m (gi,g 2 )= 9iJ n^JSX^s _ Wn«j* + rm+l^ ' ( 28 ) 



V n^+(m+l)«i[ 5,2 V n g * + (m+l) g * 

n „ / m+1 

2(n«$ + (m+l)«tf) fly 2(n<$+(m+l)<tf) 
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with corresponding Euler angles, 

0i = arccos[-H^T °2 = -arccos[ . ^ =1, 6 3 = arccos[ . £m J 



(29) 



The time- dependent probability amplitudes of the three levels are given by 



c 3 n - i,m (t) 

cr 1,m+1 (t) _ 










- iQ()t 








T n>m {gi,g2) 



c 3 n - 1,m (o) 

C 2 n ' m (0) 

cr 1,m+1 (o) 



(30) 

Now similar to the semiclassical model the probabilities corresponding to different initial 
conditions are: 



Case-IV: When the atom is initially in level-1, i.e, C 1 



n— l,m+l 



and 



C*3 1,m = 0, the time-dependent atomic populations of the three states are given by 



a 



n—l,m 



(t) 



n— l,m+l 



(0 



1 £ n 



i N| + (m + 1)# 2 cos fi„ m t] 2 . 



(31a) 
(31b) 
(31c) 



Case-V: When the atom is initially in level-2, i.e, C?- 1 ^ 1 = 0, C 2 " ' m = 1 and C 3 " _1 ' m = 0, 
the probabilities of three states are 



C 3 ' (£) — t^I 2 - sin 2 Q nm t, 



l^2 n ' m W| 2 = n*-[(™+ l)^ + n^co8fi Bm t] 2 , 



*~in— l,m+l 
U 3 



Case- VI: If the atom is initially in level-3, then we have CY 
1 and the corresponding probabilities are 

^"-^(Of^cos 2 ^, 

cr 1,m+1 (i) 



n— l,m+l 



0, c\ 



1 'nr. . 

2 _ (m+l)a2 . 2 Q . 



(32a) 
(32b) 

(32c) 
and 

(33a) 
(33b) 
(33c) 



We now proceed to evaluate the population oscillations of different levels of the vee system 
with similar initial conditions. 



11 



V.The semiclassical vee system 



At zero detuning, the Hamiltonian of the semiclassical three-level vee system interact- 
ing with two-mode classical fields is given by 



H v = 



hui h,Ki exp[— iflit] 

T1U2 h,K,2 exp[— iQ 2 t\ 

tiKi exp[iQit] hK2exp[iQ 2 t] — Ti{ui\ -\- UI2) 



(34) 



Let the solution of the Schrodinger equation corresponding to Eq.(34) is given by 

= |1) + C 2 (t) |2) + C 3 (t) |3) , (35) 

where Ci(t), C 2 (t) and C 3 (t) are the time- dependent normalized amplitudes with the 
basis vectors defined in Eqs.(13). To calculate the probability amplitudes of three states, 
substituting Eq.(35) into the Schrodinger equation we obtain 



i^jja = + Kl exp(-iftit)Ci, 



= u 2 C 2 + « 2 exp(-ifi 2 *)Ci, 



i^- — —{oj-i + co 2 )C 1 + k 2 exp(ifi 2 t)C 2 + «i exp(iQ 1 i)C 3 . 

Let the solutions of Eqs.(36) are of the following form: 

C 3 (t) = A 3 exp(iS 3 t), 

C 2 {t) = A 2 exp(iS 2 t), 
d(t) = Aiexp(iSit), 



(36a) 
(36b) 
(36c) 

(37a) 
(37b) 
(37c) 



where Aj-s are the time independent constants to be determined from the boundary 
conditions. From Eq.(36) and Eq.(37) we obtain 



(Si - fii +coi)A 3 + rexAi = 0, 
(Si - fi 2 + cj 2 )A 2 + re 2 Ai = 0, 
(Si —u 2 — u;i)A 1 + k 2 A 2 + K1A3 = 0. 



(38a) 
(38b) 
(38c) 



In deriving Eqs.(38), the time independence of the amplitudes A3, A 2 and Ai are ensured 
by invoking the conditions S 2 = Si — Q 2 and S 3 = Si — fti. At resonance, AY = = — 
i.e. (2u; 2 + u>i) — Q 2 = = (cj 2 + 2u>i) — Qi and the solutions of Eq.(38) are given by 



Si = (loi + u 2 ) 



(39a) 
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Si = {uji + u 2 )±A (39b) 

and we have three values of S 2 and £3 

S^ = -oj 2 , S 2 ' 3 = -co 2 ± A (40a) 

Si = -wi, S 2 / — -ui ± A. (40b) 

Using Eqs.(39) and (40), Eqs. (37) can be written as 

C 3 (t) = A^exp(-^it) + A 2 exp(-i(cj! + A)*) + A|(-i(wi - A)t), (41a) 

C 2 (t) = A^exp(-?cj 2 t) + A 2 exp(-?(w 2 + A)t) + A^(-?(w 2 - A)t), (41b) 

Ci(t) = exp(i(w 2 + wi)*) 

+A 2 exp(i((cj 2 + wi) + A)t) + A? exp(i((cc> 2 + wi) - A)t), (41c), 

where Aj-s are the constants which are calculated below from the various initial conditions. 

Case-I: Let us consider initially at t — 0, the atom is in level- 1, i.e, Ci(0) = 1, 
C 2 (0) = and 6*3(0) = 0. Using Eqs. (36) and (41), the time dependent probabilities of 
the three levels are given by 

|C 3 (*)| 2 = grSin 2 Ai, (42a) 

|C 2 (t)| 2 = ^sin 2 At, (42b) 

Id (t) | 2 = cos 2 At. (42c) 

Case-II: If the atom is initially in level-2, i.e, Ci(0) = 0, C 2 (0) = 1 and C 3 (0) = 0, the 
corresponding probabilities of the states are given by 

|C 3 (t)| 2 = 44#sin 4 At/2, (43a) 

\C 2 (t)\ 2 = ^(z^ + ^cosAt) 2 , (43b) 

|Ci(*)| 2 = ^ sin 2 At. (43c) 

Case-Ill: When the atom is initially in level-3, i.e, Ci(0) = 0, C 2 (0) = and C 3 (0) = 1, 
we obtain the the occupation probabilities of the three states as follows: 

|C 3 (t)| 2 = ^(z^ + ^cosAt) 2 , (44a) 
|C 2 (t)| 2 = 44#sin 4 At/2, (44b) 
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|Ci(f)| 5 



^ sin 2 At. 



(44c) 



VI. The quantized vee system 

The eigenfunction of the quantized vee system described by the Hamiltonian in Eq.(lO) 
is given by 



|*v(0) = E [Cr +1 ' m (t) \n + 1, to, 1) + C 2 n ' m (t) |n, to, 2) + C" +1 ' m -\t) \n + 1, to - 1, 3)]. 

ra,rrt=0 



Once again we note that the Hamiltonian couples the atom-field states \n + l,m, 1), 
1 72, to, 2) and |ra + 1, to — 1, 3) forming vee configuration depicted in Fig. 2. The inter- 
action part of the Hamiltonian (45) can also be expressed in the matrix form 



Hg = h 



9W m 




g 2 Vn + T 

giVm g 2 Vn + 1 



(46) 



and the corresponding eigenvalues are A± = zLh^rngf + (n + \)g\ (= ±^nm) and Ao 
respectively. The dressed eigenstate is given by 



\nm, 


3) " 




\n + 1, TO — 1, 3) 


\nm, 


2) 


— T 


71, TO, 2) 


\nm, 


1) . 




71 + 1, TO, 1) 



the rotation matrix is found to be 



T = 

± n,m 



gi 



2((n+l)g%+mg'j) 



2{(n+l)gl+mgl) 



n+1 
i+l)gl+rng{ 



n ! ~~ -L 

i/2 y 2((n+l)<£ +mpjf) >/2 

^ {n+i)7j+^gj 

_n ./ n+1 ~ J_ 

y2 v / 2 ((„ + i)^ +m ^) ^ 



The straightforward evaluation yields the various Euler angles are 



0i = 



4' 



# 2 = arccos[- 



/ra+lff2 



1,^3 = 



\Jmgl+{l+n)g\ 

The time-dependent probability amplitudes of the three levels are given by 

D —iCl nm t n Q 



c n+l,m-l^ 

c 2 n - m (0 

C? +1>m (*) 



-l 










-,—iQ.ot 









T 



C 2 n ' m (0) 

cr +i ' m (o) 



(47) 



(48) 



(49) 



(50) 
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Once again we proceed to calculate the probabilities for different initial conditions. 

Case-IV: Here we consider initially the atom is in level-1 i.e, (7™ +1>m — i ; (™ > m — o 
and Cz +1,m 1 = 0. Using Eqs.(49) and (50) the time- dependent probabilities of the three 
levels are given by 



.^n+l,m— 1 
°3 



(t) 



Sill i 'nm^; 



|C*2 ,m (t) | — - i ^ 2 sin 2 f2 nm t, 



cr +i - m (t)" = cos 2 fi nm t. 



(51a) 
(51b) 
(51c) 



Case-V: If the atom is initially in level-2 i.e, C£ +1 ' m_1 = 0, C 2 n,m = 1 and C{ l+1 ' m = 0, 
then corresponding probabilities are 



s-yn+l,m— 1 , 
°3 



(tr = 4 ^(n+l)M sin 4^ mnt/2; 

™ L m.n. 



\CZ m {t)\ 2 = ^[mgj + (n + l)^cosfi„ m t] 2 , 



(52a) 
(52b) 

(52c) 



Case- VI: Finally if the atom is initially in level-3 i.e, C? +1 ' m = 0, C 2 n ' m = and C 3 n+1 ' m " 1 = 
1, then 



\cr n {t)\ 2 = 4 ^?(yi)M ^n nm t/2, 



C\ + ' (t) — qt 1 - sin 2 VL nm t. 



(53a) 
(53b) 
(53c) 



Finally we note that for large values of n and m, Case-IV, V and VI become identical to 
Case-I, II and III, respectively. This precisely shows the validity of the Bohr's correspon- 
dence principle indicating the consistency of our approach. 



VII. Numerical results and discussion 

Before going to show the numerical plots of the semiclassical and quantized lambda 
and vee systems, we first consider their analytical results. If we compare Case-I, II, III 
of both cases, we find that the probabilities in Case-I (Case-Ill)) of lambda system is the 
same as in Case-Ill (Case-I) of vee system except the populations of 1st and 3rd levels 
are interchanged. See Eqs.(21 & 44) and Eqs.(23 & 42) for detailed comparison. Also 
Case-II respective models are similar which is evident by comparing Eqs.(22 & 43). In 
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contrast, for the quantized model, Case-IV (Case- VI) of the lambda system is no longer 
same as in Case- VI (Case-IV) of the vee system. This breaking of symmetry is evident 
by comparing the analytical results, Eqs.(31 & 53), Eqs.(32 & 52) and Eqs.(32 & 51) 
respectively. Unlike previous case, also Case-V both the models are different which is 
evident from Eqs.(22 & 43). 

In what follows, we compare the probabilities of the semiclassical and quantized 
lambda and vee systems respectively. Fig. 3 and 4 show the plots of the probabilities 
|C{(t)| 2 (blue line), \C % 2 (t)\ 2 (green line) and | C3 (^) | 2 (red line) for the semiclassical lambda 
and vee models when the atom is initially at level-1 (Case-I), level-2 (Case-II) and level-3 
(Case-Ill), respectively. The comparison of the plots shows that the pattern of the prob- 
ability oscillation of the lambda system for Case-I shown in Fig.3a (Case-Ill in Fig.3c) 
is similar to that of Case-Ill shown in Fig.4c (Case-I in Fig.4a) of the vee system. More 
particularly we note that in all cases the oscillation of level-2 remains unchanged, while 
the oscillation of level-3 (level-1) of the lambda system for Case-I is identical to that of 
level-1 (level-3) of the vee system for Case-Ill. Furthermore, comparison of Fig.3b and 
Fig. 4b for Case-II shows that the time evolution of the probabilities of level-2 of both 
systems also remains similar while those of level-3 and level-1 are interchanged. From 
the behaviour of the probability curve we can conclude that the lambda and vee config- 
urations are essentially identical to each other as we can obtain one configuration from 
another simply by the inversion followed by the interchange of probabilities. 

For the quantized field, we first consider the time evolution of the probabilities taking 
the field is in a number state representation. In the number state representation, the 
vacuum Rabi oscillation corresponding to Case-IV, V and VI of the lambda and vee 
systems are shown in Fig.5 and Fig. 6 respectively. We note that, unlike previous case, the 
Rabi oscillation for Case-IV shown in Fig. 5a (Case- VI shown in Fig.5c) for the lambda 
model is no longer similar to Case- VI shown in Fig.6c (Case-IV shown in Fig. 6a) for the 
vee model. Furthermore, we note that for Case-V, the oscillation patterns of Fig.5b is 
completely different from that of Fig.6b. In a word, for the quantized field, in contrast 
to the semiclassical case, the symmetry in the pattern of the vacuum Rabi oscillation in 
all cases is completely spoiled irrespective of the fact whether the system stays initially 
in any one of the three levels. 

The quantum origin of the breaking of the symmetric pattern of the Rabi oscillation 
is the following. We note that due to the appearance of the terms like (n + 1) or (m + 1), 
several elements in the probabilities given by Eqs. (31,32,33) for the lambda system and 
Eqs. (51,52,53) for the vee are non zero even at m = and n = 0. We argue that the 
vacuum Rabi oscillation interferes with the probability oscillations of various levels and 
spoils their symmetric structure. Thus as a consequence of the vacuum fluctuation, the 
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symmetry of probability amplitudes of the dressed states of both models formed by the 
coherent superposition of the bare states is also lost. In the other word, the invertibility 
between the lambda and vee models exhibited for the classical field disappears as the 
direct consequence of the quantization of the cavity modes. 

Finally we consider the lambda and vee models interacting with the bi-chromatic 
quantized fields which are in the coherent state. The coherently averaged probabilities of 
level-3, level-2 and level-1 are given by 

2 



{P 3 (t)) A = E W n W m 

n,m 

(P2(t)) A = E W n W n , 

n,m 

(Pi(t)) A = E W n W m 

n,m 

for the lambda system and 

(P 3 (t)) v = E W n W„ 

n,m 

{P*(t)) v = E W n W n 

n,m 

(Pi(t)) v = E W n W n 



n—l,m 



a 



I (j n - ) 



(!) 



(Ol 2 , 



n— l,m+l 



(t) 



fyn+l,m— 1 



(t) 



\^jri,m 



(t)\ 2 , 



C[ l+1 > m (t) 



(54a) 
(54b) 

(54c) 

(55a) 
(55b) 

(55c) 



for the vee system, where W n = ^,exp[— n]n n and W m = ^jexp[— m]m m with n and fh 
be the mean photon numbers of the two quantized modes, respectively. Fig. 7-9 display 
the numerical plots of Eq.(54) and Eq.(55) for Case- IV, V and VI respectively where the 
collapse and revival of the Rabi oscillation is clearly evident for large average photon 
numbers in both the fields. We note that in all cases, the collapse and revival of level-2 of 
both the systems are identical to each other. Furthermore, we note that the collapse and 
revival for lambda system initially in level-1 shown in Fig.7a, Fig. 7b and Fig. 7c (level-3 
shown in Fig. 9a, Fig. 9b and Fig. 9c) is the same as that of the vee system if it is initially 
in level-3 shown in Fig.7f, Fig.7e and Fig.7d (level-1 shown in Fig.9f, Fig.9e and Fig.9d) 
respectively. On the other hand, if the system is initially in level-2, the collapse and 
revival of the lambda systems shown in Fig.8a, Fig. 8b and Fig.8c are identical to Fig.8f, 
Fig.8e and Fig.8d respectively for the vee system. This is precisely the situation what 
we obtained in case of the semiclassical model. Thus the symmetry broken in the case 
of the quantized model is restored back again indicating that the coherent state with 
large average photon number is very close to the classical state where the effect of field 
population in the vacuum level is almost zero. It is needless to say that the coherent state 
with very low average photon number in the field modes can not show the symmetric 
dynamics in lambda and vee systems. 



17 



VIII. Conclusion 



This paper presents the explicit construction of the Hamiltonians of the lambda, vee 
and cascade type of three-level configurations from the Gell-Mann matrices of SU(3) group 
and compares the exact solutions of the first two models with different initial conditions. 
It is shown that the Hamiltonians of different configurations of the three-level systems 
are different. We emphasize that there is a conceptual difference between our treatment 
and the existing approach by Hioe and Eberly [18,21,22]. These authors advocate the 
existence of different energy conditions which effectively leads to same cascade Hamilto- 
nian (/i 2 i 7^ 0, h 32 ^ and /i 31 = in Eq.(l)) having similar spectral feature irrespective 
of the configuration. We justify our approach by noting the fact that the two-photon 
condition and the equal detuning condition is a natural outcome of our analysis. For 
the lambda and vee models, the transition probabilities of the three levels for different 
initial conditions are calculated while taking the atom interacting with the bi-chromatic 
classical and quantized field respectively. It is shown that due to the vacuum fluctuation, 
the inversion symmetry exhibited by the semiclassical models is completely destroyed. 
In other words, the dynamics for the semiclassical lambda system can be completely ob- 
tained from the knowledge of the vee system and vice versa while such recovery is not 
possible if the field modes are quantized. The symmetry is restored again when the field 
modes are in the coherent state with large average photon number. Such breaking of the 
symmetric pattern of the quantum Rabi oscillation is not observed in case of the two-level 
Jaynes-Cummings model and therefore it is essentially a nontrivial feature of the multi- 
level systems which is manifested if the number of levels exceeds two. This investigation 
is a part of our sequel studies of the symmetry breaking effect for the equidistant cascade 
three-level and equidistant cascade four-level systems respectively [36,37]. Following the 
scheme of constructing of the model Hamiltonians, it is easy to show that we have different 
eight dimensional Bloch equations and non-linear constants for different configurations of 
the three-level systems and these issues will be considered elsewhere [40]. The breaking 
of the inversion symmetry of the lambda and vee models as a direct effect of the field 
quantization is an intricate issue especially in context with future cavity experiments with 
the multilevel systems. 
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Fig. 2 : Vee type transition 



22 



(a) (b) (c) 




t t t 

[Fig. 3]: The time evolution of the probabilities |Ci(t)| 2 (blue line), | C2 (*) 1 2 (green 
line) and | C3 (i ) | 2 (red line) of the semiclassical lambda system for Case-I, II and III 
respectively with values k\ = .2, K2 = •!■ 



(a) (b) (c) 




t t t 

[Fig.4]: The time variation of the probabilities |Ci(£)| 2 (blue line), | C2 (^) 1 2 (green 
line) and | C3 ) | 2 (red line) of the semiclassical vee system for Case-I, II and III 
respectively with above values of K\, k^. 
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(a) (b) (c) 




t t t 

[Fig.5]: The Rabi oscillation of the quantized lambda system when the fields are 
in the number state for Case-I, II and III, respectively with gi = .2, g2 = .1, n = 1, 
m = 1. 



(a) (b) (c) 




t t t 

[Fig. 6]: The Rabi oscillation of the quantized vee system when the fields are in the 
number states for Case-I, II and III, respectively for same values of gi, g2, n, m. 
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[Fig. 7]: Figs.7a-c display the time-dependent collapse and revival phenomenon of 
level-3, level-2 and level-1 of the lambda system for Case-IV, while Figs.7d-f show 
that of the level-3, level-2 and level-1 respectively for of Case- VI of the vee system 
taking the field modes are in coherent states with n = 30 and m = 20. 
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[Fig. 8]: Figs.8a-c display the time-dependent of collapse and revival of level-3, level- 
2 and level-1 of the lambda system for Case-V while Figs.8d-f show that of level-3, 
level-2 and level-1 of the vee system for Case-V with the same values of n and fa. 
as in Fig.7 
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[Fig. 9]: Figs.9a-c display the time-dependent of collapse and revival of level-3, level- 
2 and level-1 of the lambda system for Case- VI while Figs.9d-f show that for level-3, 
level-2 and level-1 respectively for the vee system for Case-IV with the same values 
of n and fa as in Fig. 7. 
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